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Abstract 

We present a mathematically rigorous method suitable for solving three-body 
bound state and scattering problems when the inter-particle interaction is of 
a hard-core nature. The proposed method is a variant of the Boundary Con- 
dition Model and it has been employed to calculate the binding energies for 
a system consisting of three ^He atoms. Two realistic He-He interactions of 
Aziz and collaborators, have been used for this purpose. The results obtained 
compare favorably with those previously obtained by other methods. We 
further used the model to calculate, for the first time, the ultra-low energy 
scattering phase shifts. This study revealed that our method is ideally suited 
for three-body molecular calculations where the practically hard-core of the 
inter-atomic potential gives rise to strong numerical inaccuracies that make 
calculations for these molecules cumbersome. 
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I. INTRODUCTION 



The Boundary Condition Model (BCM) (see, for example, Refs. [0,0]) is of interest 
due to its simplicity in describing the short-range component of the interaction between 
particles. In the BCM the interaction is specified by boundary conditions imposed on the 
wave function when the particles approach each other at a certain distance r = c. The so- 
called hard core potentials represent a particular variant of the BCM where one requires that 
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the wave function vanishes at r = c. Such a requirement is equivalent to an introduction of 
an infinitely strong repulsion between particles at distances r < c. The standard formalism 
for scattering [|,|| does not deal with hard-core interactions described by these boundary 
conditions. Therefore a derivation of special equations to handle this class of interactions is 
desirable. 

Replacement of the finite, for r > 0, but often singular at r = 0, repulsive short-range 
part of the potential with a hard-core interaction turns out to be a very effective way to 
suppress inaccuracies related to a numerical approximation of the Schrodinger operator at 
short distances. Although in two-body applications these potentials are easy to handle, 
in three-body systems certain mathematical difficulties appear 0, which are absent when 
conventional potentials are used. 

To overcome these difficulties various approaches were considered. We shall recall here 
the two main ones related to the Faddeev equations 0,^. In the first, a certain limiting 
procedure is used where special potentials that include only a finite repulsive core are con- 
structed at a first step. The parameters of these potentials are then chosen so that the final 
two-body wave function satisfies the desired boundary conditions ||^-§|. The corresponding 
two-body t-matrices are subsequently substituted into the Faddeev integral equations 
under an implicit assumption that the latter are still valid. The resulting equations are con- 
sidered as a generalization of the Faddeev equations for the BCM. A similar approach was 
also used in Refs. ||I1,|9|JT0| . A common feature of the reduced equations in these approaches is 
that they are not of a Fredholm type and that they have not a unique solution at all energy 
values, including the complex ones. To obtain a unique solution one is forced to introduce 
auxiliary conditions or relations p|,p!0|. 

In the second approach, three-body integral equations of a Fredholm type are derived 
in the BCM model without any limiting procedure. Instead, one uses the fact that the 
spectral problem for the Schrodinger operator is an example of a classical boundary- value 
problem for elliptic differential equation in partial derivatives of the second order. One of the 
traditional methods to deal with such problems is the Potential Theory [|1^ . An approach to 
the three-body problem in the BCM which is based on the Potential Theory was developed 
in Refs. |T2HI5[ (see also and [0). However, in contrast to the boundary value problems 
for compact surfaces, the initial three-body equations were not of Fredholm type, similarly 
to the three-body Lippmann-Schwinger equation in the case of the conventional soft-core 
potentials. This is due to the same reason, i.e, the noncompactness of the support of the two- 
body interaction in the three-body configuration space. To overcome this problem, these 
equations were rearranged in Refs. [p!2|-[T5|| using the Faddeev method IQ. The resultant 
equations are of Fredholm type and suitable for use in the three-particle scattering problem. 
These developments allowed the reformulation of the Faddeev equations for the bound state 
and scattering problems in configuration space in terms of boundary-value problems which 
are suitable for numerical calculations. This was demonstrated in three-nucleon bound-state 
and scattering calculations |]13| , |16| . 



In this work we shall present a hard core version of the BCM formalism [P^ T3] and apply 
it to the three-atomic "^He system. Various methods have been used in the past to study the 
ground state properties of '^He molecules. We mention here the Variational Method (VM) 
T^-gOl, the Variational Monte Carlo method (VMC) the Green Function Monte 



Carlo method (GFMC) [P2|-|27|, the methods based on the Faddeev integral equations in 
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momentum space p8|j29[] , the Faddeev differential equations in configuration space [Q, and 



the hyperspherical methods ||3T|-p3 



The general atom-diatom collision problem has been addressed by various researchers 
in the field and we refer the interested reader to the review articles on this topic by Micha 
3^ and Kuppermann |^5[. Collision dynamics at thermal energies of the H+H2 system and 



the existence of resonances were discussed by Kuruoglu and Micha |^ using the Faddeev 
integral equations in momentum space. Finally, the problem of existence of the ^He n-mers 
and its relation to the Bose-Einstein condensation in He II was discussed in Refs. []37| , p8 . 
From the experimental works we mention those of Refs. where molecular clusters 



consisting of a small number of noble gas atoms were investigated. 

The interaction between bosons in such clusters is usually described by central potentials 
having a very strong repulsive cores p3|-|47[|. In the present work, we approximate the 
strong repulsion between the Helium atoms at short distances by a hard core and solve the 
corresponding boundary value problems for the Faddeev-type differential equations. We 
shall show that the method gives excellent results for the ground-state energy of the Helium 
^He trimer. It further allowed the calculation of an excited state interpreted in p9|j33[1 as an 
Efimov one [|^. Moreover, we shall demonstrate that the method is suitable for scattering 
calculations at ultra-low energies below as well as above the breakup threshold. Certain 



results of our work were presented in |H9 



Some comments on the notation used throughout the paper: The y^, z 
for the main branch of the function z^l"^ . The a denotes the unit vector, a = 



G C, stands 

A, a G M", 
|a| 

while L2{D) is the standard notation used for the Hilbert space of the square integrable 
functions defined in a domain D of M". The symbol W2{D) stands for the space of those 
of the functions of L2{D) which have all second partial derivatives as elements of L2{D). 
Finally, the notation D is used for closure of a set D C M". 

This paper is organized as follows. In Sec. II we overview the three-body bound and 
scattering state formalism for the hard core interactions and in Sec. HI, we describe its 
application to a system of three identical bosons. Our numerical results for the three-atomic 
^He system are presented in Sect. IV while our conclusions are drawn in Sec. V. A detailed 
description of the numerical methods used is given in the Appendix. 



II. THREE-PARTICLE SYSTEMS WITH HARD CORE INTERACTIONS 



a 



In describing the three-body system we use the standard Jacobi coordinates Xq,, ja, 
= 1,2,3, expressed in terms of the position vectors of the particles G and their 



2m,m, 1 
mp + m^_ 

2mQ,(m^ + m^) 
m^ + m/3 + m^_ 



1/2 



(1) 



r„ - 



m/3 + m^ 



where (a, 7) stands for a cyclic permutation of the indices (1, 2, 3). The coordinates x^,, y„ 
fix the six-dimensional vector X = {^ajyaj £ The vectors x^, corresponding to the 
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same point X as the pair ^a^Ya are obtained using the transformations 

where the coefficients c^a and s^^ fulffil the conditions — 1 < c^q, < +1 and s^^ = 1 — c^^ 
with Ca/3 = C/3q, Saf3 = —SjSa, jS ^ a and depend only on the particle masses Q. For equal 
masses C/jq = —1/2. 

The configuration space Q of the three-body system in the hard-core model represents 
only a part of the space MP external, |xq,| > c^, with respect to all three cylinders Fq,, 
Fq, = {X G : X = (xa,ya), |xq,| = Ca}, a = 1,2,3, where Cq, > 0, stands for the 
values of |xq| when the cores of the particles in the pair a contact each other. A three- 
dimensional image of this space for particles with equal masses and the same core radii 
Cq, = c, a = 1, 2, 3, c > 0, is sketched in Fig. 1, in coordinates Xa = |xq,|, i/a = \ya\, and 
rja = ^aYa- The cylinders Fi, F2 and F3 are depicted in this figure by the plane Xi = c and 
surfaces jxl + jyf—^xiyirji = c^, jxf + jyl+^xiyirii = c^, respectively. The domain f2 is 
a part of the set > c, ?/i > 0, — 1 < ?7i < 1, restricted by F2 and F3. From a geometrical 
point of view, the image shown in Fig. 1 is exactQ since only coordinates (such as Eulerian 
angles) describing a rotation of the plane defined by the position of particles are omitted. 

The Hamiltonian of a system of three particles with hard-core interactions is defined in 
L2(f2) by the expression 

Hf{X)= (^-^x + f2v}jf{X) (2) 

on the set of functions f{X), f G 1^2^ (fi), satisfying the condition 

/ Ian = (3) 

on the boundary dQ of the domain Q. The Laplacian —Ax corresponds to the kinetic energy 
operator of the system under consideration. The potentials Va, a = 1,2,3 are two-body 
interactions and thus when acting on the function / in the expression (^ they only operate 
on the corresponding two-body variable x^, |xq| > c^. We assume that these pair potentials 
are bounded Hermitian operators. The Hamiltonian H that includes such potentials is a 
self-adjoint operator and thus its spectrum is real. For local potentials we assume that 

|K.(x„)| <Cq(1 + |x„|)-3-% |x„| >c„, (4) 

where the constants Cq > and e > 0. Similar conditions are assumed for the partial deriva- 
tives of V"q,(xq,). The Aziz et al. potentials considered in this work are examples of 



such interactions. Nonlocal potentials can also be included in our formalism, provided their 
kernels Vq(xq,, x'q) are smooth functions obeying conditions similar to (H) 



It should be noted that the transition to the variables Xq,, ya and t/q, is not conformal. In 
particular the true angle between any two surfaces Tp and F^, /3 7^ 7, at points belonging to the 
intersection manifold F^PlP^y varies between ^ — (fi^^ and ^ + (/>^^, 4)^^ = arcsin|c^^| and never 
acquires the value of or vr (since |c/3^| < I). 
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A. Bound state problem 



We shall consider first the boundary value problem for the Faddeev differential equations 
for the three-body bound state. Let 

i/^ = , 

E being the bound state energy and \Ef the corresponding three-body bound state wave 
function. We are concerned with states for which E < and that these energies are below 
the threshold of the continuous spectrum of H. Using the Green's formula (see, e.g., Ref. 
|ll[|) one can easily show that the function satisfies the following Lippmann-Schwinger 



type equation 

^{X) = - j dosG,{X,S-E)^^>{S)-Y, j dX'Go{X,X'-E){V^m){X') (5) 
an n 

where Gq{X,X'\ z), is the three-body free Green function, i.e., the kernel of the resolvent 

G,{z) = (-Ax - z)-^ 

of the Laplacian —Ax in the six-dimensional space M^. We recall that the function 
Gq{X^ X'; z) can be expressed in terms of the Hankel function of the first kind ifj^^ 

The ns denotes the external unit vector (directed into Vt) normal to the surface dVL while 
das is a surface element (five-dimensional square) on dVt. 

In contrast to '^{X), defined only for X G fi, the function Go{X, X'\ E) is defined for all 
X G M^, X 7^ X'. Therefore, the right-hand side of (|^) is defined for X G as well as for 
X G \ fi. Moreover, from the Green's formula it follows that, for any X G \ 



j dasGo{X,S;E)^^{S)-J2 J dX' Go{X, X'; E) {V^^){X') = . 

Qr» ^ a=l r-t 



(6) 



an 

The Faddeev components of the function \E' are introduced via the formulas (see Refs. 



mm) 

$,(X) = - J dasG,{X,S-E)^^>{S)- j dX' G,{X,X'-E) {V^m){X') . (7) 



We shall consider the functions $a(X) given by (^ for all X G M®, i. e., outside as well as 
inside the surface dVL. From (Bl) and (||) one gets 



2^<i>„(x)-| XGR6\n. 

a=l ^ ^ 



(8) 
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The surface integral 



dasGo{X,S;z) fia{S) , zE 



(9) 



which appears in (0), represents the potential of a simple layer |]TT| with density /Iq, concen- 
trated on the surface Fq,. In our case 



d 

dns 





lis eV^\dVl. 



As has been shown in Refs. p2|-p!^, each of the densities fia{S) on the cylinder as a 
function of the variable S" G Fq,, is everywhere continuous except perhaps where this cylinder 
intersects the other two cylinders F^, (3 ^ a. This means that in crossing the surface Fq (at 
least not on the intersection of Fq with F/3), the potential of a simple layer (|^) is a continuous 
function Evidently for X ^ Fq, the integral is infinitely differentiable with respect 
to X and that 

(-Ax -z) f das Go{X, S; z) fi^{S) = . 



Acting on both sides of the equality (|^ by the differential expression —Ax — E and taking 
into account the relation (H), one obtains the following system of differential equations for 
the components $q,(X), 



-Ax + V^-E)<^^{X) = -V^Y.<i>fs{X), |xj>c. 



(10) 



-Ax-E)^^{X) = 



X/-y ^ Cfy . 



According to (H), the sum of the functions $a(X) must vanish not only on the surface dQ 
but also inside of it, i.e., 



^$^(x) = o, XeR^\n. 
0=1 



(11) 



In fact one can replace the very strong conditions (pJ]) with the essentially more weak 
conditions |1I2| . |T3[ 



/3=1 



a = 1, 2, 3 , 



(12) 



l^a I — Cq 



requiring that the sum of $q,(X) to be zero only on the cylinders Fq,. It is understood that 
for the bound-state problem, the conditions 



G L^iR"), 



a = 1, 2, 3 , 



(13) 
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must be fulfilled. Similarly to the pure potential model the asymptotic behaviour of 
as |X| — i> oo is of an exponential character, the form of which is quite complicated 0. 

Equations (p!OD, (]I2|), and (|13|) describe the boundary value problem for three-body bound 
systems with hard-core interactions and are a natural generalization of the Faddeev differ- 
ential formulation [^] for bound states. 

The numerical advantage of our approach is already obvious from the structure of Eqs. 
(PUj): When a potential with a strong repulsive core is replaced with the hard-core model, one 
approximates, inside the core domains, only the Laplacian — A^, instead of the sum of the 
Laplacian and a huge repulsive term, and in this way a much better numerical approximation 
can be achieved. 



B. Scattering processes 

Let ^['3'«l±(X,p^) be the three-body wave function corresponding to a (2 -|- 1 ^ 2 -|- 
1 ; 1 -|- 1 -|- 1) process where in the initial state the pair subsystem (3 is bound in a state 
with energy e^j^^, e^^g < 0, and the complementary particle is asymptotically free, the 
relative momentum being p^, G M'^. By ^ we denote here a distinctive label (consisting of 
appropriate quantum numbers) for the two-body state concerned. The Faddeev components 
|f|,|l6| $a(X) = <l>[f'^'^(X,p^) of the wave function ^[^'«]±(X, p^), 

vl/[A€l±(X) = ^$[^'«l±(X), 

a=l 

in the hard-core model satisfy the same differential equations ([ToD and boundary condi- 
tions ([T^ of the three-body bound state problem. These components can be written as 

^f^^Hx, vp) = 5a^x,Ax, Pp) + Yl M^^)ul':tiya, pp) + ulf^ix, p^) (i4) 

where 



XpAX, P/?) = V^/3,5(x/3) exp(ip^-y^) 



(15) 



is the incident wave consisting of a two-body bound state tp^^^ and a plane wave. The 
functions uj^'^^'^ and uj^Q^"^ have the same asymptotic behavior [|15[ as in the usual potential 
model M, namely. 



±iy^E-e 



"|y<.| 



U, 



[/3,C]± 



a,0 



ya ^ oo 



X oo 



,±iv^|X| 
|X|5/2 



aL^|/'^(ya,P/3) + o(|yc 



4^'^^^(^,p,)+o(|xrV2) 



-1/2N 



(16) 
(17) 



where E = e^^^ + P% is the energy of the system. For E > ea^^i the function a^^^,'^(yQ, p/3) 
represents the amplitude for the elastic (a = /3, ^' = ^) or rearrangement {a ^ (3 oi E,' ^ 
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scattering. The functions '^'^(X, p^) provide us, at -E > 0, with the Faddeev components 
of the total breakup amphtude ^[^'^'^(X, p/j) 

It should be stressed that the two-body eigenfunctions ipa,^{'^a) are assumed to be zero 
within the respective cores a, i.e., ipa,c{^a) = for |xa| < c^. The boundary- value problem 
as described by Eqs. (|10|), (0), and (|14|)-(0), is the extension of the Faddeev formalism to 
the (2 + 1^2-1-1; 1-|-1-|-1) scattering processes for hard-core potentials. 

A detail analysis for the boundary-value problems described above, the derivation of 
the asymptotic boundary conditions for scattering states as well as other boundary-value 
formulations, can be found in Refs. Here, we only recall, briefly, some peculiar 

properties of the discrete spectrum generated by the condition (|1^). As compared to the 
spectrum of the initial Hamiltonian defined by Eqs. (Q) and (^, this spectrum acquires 
an additional component, corresponding to the Dirichlet boundary- value problems for the 
domains which result from the intersection of the cylinders Fq. We introduce the following 
notations for these domains: Let Aq,/?^ be a domain restricted by all the three cylinders Fq, 
a = 1, 2, 3, and (9Aq^^ be its boundary (see Fig. 1). The notation Aa/s is used for a part of 
the domain bounded by the cylinders Fq, and F/3, (3 7^ a, and at the same time is external 
with respect to the set Aq,^^. By Aq, we denote the domain bounded by the cylinder F^ 
which is at the same time external to the rest cylinders F/3, (3 ^ a. The notations dhap and 
dha are used for the boundaries of the domains K^p and A^, respectively. 

It can be shown [|T^] that the discrete spectrum of the boundary- value problem (jl^, 
and ( |T^ ) includes not only the discrete spectrum cXdiH) of the original Hamiltonian H but 
also a set cij^"^ consisting of a discrete set of eigenvalues of the homogeneous internal Dirichlet 
problems in the domains A^/j^, Aap, and A^, a,P = 1,2,3, /? 7^ a, namely, the discrete 
spectra of the operators defined in W2{Aai3'y), W2{Aa), and W2{Aai3) by the expression (H) 
(under the assumption (Vq,/)(xq,) = for |xq| < Cq) and the respective boundary conditions 
f \dA ,3 ~ 0; / I^A^ ~ f \dA ,3 ~ ^- There exists a simple criterion in selecting solutions 
of the spectral problem described by (p!0D , (IT^), and ([T3|) corresponding to the spectrum of 
the Hamiltonian H only. This is just the requirement (see also the condition (|TT])) that the 
total wave function inside the cylinders F^ vanishes, 

3 

\If{X) = ^^fs{X) = if|x„|<c^, a = 1,2, 3. 

13=1 

It should be noted that the lower boundary inf a'^'^^ of the auxiliary spectrum a^"^ is situated 
above0 the lower boundary of the spectrum of the Hamiltonian H. Therefore in searching 
for a ground state no validity check of this criterion is necessary. 

The elements of the set cr^*^^ are points where the (2 + 1— i>2 + l;l + l + l) scatter- 
ing problems (p!oD , (|T^) , (p!^)-([T7|) have no a unique solution [|15|. However the auxiliary 



^In the case of sufficiently small Ca, inf crj^"^ is positive and behaves as where c = max Cq 



a 
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spectrum is discrete and thus in practice a coincidence of the scattering energy E with 
a point of the set cr^"^ can be considered as an exceptional case. In principle there is a 
way to avoid such a coincidence namely by shifting the spectrum cr^"^. This can be made, 
for example, by replacing the zero values of the potentials Va inside the core domains by 
appropriate positive values. Such a replacement does not affect the total wave function 
= in the physical domain, that is, at X G fi. 



III. PARTIAL BOUNDARY- VALUE PROBLEMS 

In what follows we shall concentrate on a system of three identical bosons interacting via 
a central potential V , i.e., via V^(xq,) = V(|xq,|), a = 1,2, 3. The total wave function of the 
system is invariant under the permutation of particles belonging to any pair a, P^^l/ = 
where is the permutation operator. This means that XcYq.) = \l/(Xfj, Yq,), a = 
1, 2, 3. Thus from the definition of the Faddeev components (|^) one obtains 

Pa^a = (18) 

i. e. 

$a(-Xa, ya) = *a(Xa, Ya) ■ (19) 

Furthermore 

$^ = P+<I>,, <I'-, = P-$„ (20) 

where P^ stand for operators of cyclic permutation of particles 

P+(123) = (312), P-(123) = (231). (21) 

The conditions (^) mean that the total wave function "^{X) is written as 

^ = (/ + P+ + p-)<l'a (22) 

where / is the identity operator. Similarly, the Faddeev equations ([T0|) and the hard-core 
boundary conditions (|T^) are written as 

{-Ax + Va- E)^^{X) = -K(P+ + P-)$a(X) , |x«| > c„ , (23) 

{-Ax-E)^^{X) = 0, |x„|<c,, (24) 

and 

<I>,(X) = -(P+ + P")<I>,(X), |xj=c„, (25) 

where ci = C2 = C3 = c and, say, a = 1. In what follows we shall drop, for convenience, the 
identification a. If one searches for a bound state of the system, the condition 

$ G L2(R^) (26) 
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is required. 

Consider now a (2 + 1 — > 2 + 1; 1 + 1 + 1) scattering process for the three bosons in an 
initial state 

X^(X,p) =^^(x)exp(ip-y). (27) 

Since the particles are identical, the incident wave X£.{^y p) must be included now, in contrast 
to ([T^), into all three summands of the total scattering wave function ^'(X) = \E'^^(X, p) 
given by equation (^) with $(X) = Therefore the Faddeev components have the 
form 

$«=^(X,p) = xdX,p) + J]^e(x)f/f (y,p) + UtiX,p) (28) 
where the terms U^^ and Uq^ have the same asymptotic form as (|T6|) and (p!7|), 

4^(y,P) = p [af (y,p) + o(|yri/2)] , (29) 

y ^ oo ^ 

±iVE\X\ I- 

^o^^(^,P) = -TJa^[^^H^^P)+o{\X\-'/')] (30) 

X oo ' ' 

where E = -\- |pp. If > e^', the function a|^(y, p) represents the elastic scattering 
amplitude, ^' = ^, or the rearrangement one, C,' ^ C,- At E > the function represents 
the Faddeev component of the total breakup amplitude ^^^(X,p) which is expressed via 

A^^iX, p) = (/ + P+ + P-) A«±(X, p) . 

The description for the auxiliary spectrum of the boundary- value problems (p3D -(p5D, 
(p6|) and (p^)-(|25|), (p8D-(pO|) is the same as the one outlined in Sect. II except that all core 
sizes are now equal, Cq = c, a = 1, 2, 3. Further we consider the case of the ^f^"*" scattering 
wave functions and thus the index "+" will be omitted. 

Similarly to Eqs. ([T0|) and (0), Eqs. (p^)-(|25|) are six-dimensional. Therefore we may 
use, for their partial wave expansion, the bispherical basis 

\IXL)= (/mA/ilLM) rr(x)lT(y) (31) 

m+/i=M 

where L is the total angular momentum of the system, Yj'^(x) and Yj^{'y), are the spherical 
harmonics, and {lmXfi\LM) the Clebsch-Gordan coefficients. 

The potential V, being central, is diagonal in the basis ( plD and has the same diagonal 
elements in all partial waves. Since the operator of the total angular momentum L and its 
projection commute with both the Laplacian —Ax and the sum I + + P~ , the study 



of the boundary- value problems (p3D -(P5D, (pBD , and (|23|) -(P5D, (P8D-(pO|) is reduced to a 
study in subspaces corresponding to fixed values of the momentum L and its projection M. 
Since the index M does not effect the structure of the equations it will be omitted. Thus 
$l(X) denotes the partial components of $(X). 
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Expanding the function $l(X) in a series of bispherical harmonics 

$^(X) = ^ ^^^^^^^ |aL) , a = {l,X}, x = |x| and = |y| , (32) 

and using the results of Ref. |Q (see also @,[16|) one obtains for (p^)-(p^) the following 
partial equations 

(ff.- = {-'''"'*;'<"'^'-^;^ (33) 

where 

^ Ijl + l) ^ XjX + l) _ 

Qrj.2 Qy2 j,2 y2 

The function "^aiix, y) represents the partial component of the total wave function ( p2D and 
is related to the partial Faddeev components ^aL^x, y) by 



^aL{x,y) = <^aL{x,y) + ^ j df] h^^, {x , y , T]) (l>a' l{x' , y') 



(34) 



where 



/ / 1 2 3 2 



,89 1 2 

y = 1/ 4 ^ + 4 ^ +^xyr], 



with T] = x-y. The functions h^^, are given by (see also 0) 

C = (-1)'^^ ^'^^,\^|f^'^ [(2A)!(2/)!(2A' + l)(2r + 1)]^/^ 
X y z 

X 5^(-i)^2A:+i)p,(.^) y: ^^7v^(-iyHv^)'^^" 

fc=0 Al+A2=A, ^ 

li+h=l 

X [(2A0!(2/i)!(2A2)!(2/,)!]-^/^ 5^(2A" + 1)(2/" + 1) ( |; ) (35) 

X"l" ^ ^ 

Ai A2 A 



A2 k l"\ k X" X' \ ( k I" r \ j I' X' L 
j I j I j I A" /" A; 



h k I 
X" I" L 



_l + X + l' + X' 



2 

where Pk{v) is the Legendre polynomial of order k. In the above, the standard notation for 
the 3-j, 6-j, and 9-j Wigner symbols, as defined in is used. It should be noted that the 
kernels h^^, depend only on the hyperangles 
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= arctan — and 9' = arctan — 
X x' 



(36) 



and not on the hyperradius 



P 



(37) 



Due to ([TsD only the components corresponding to a = {/, A} with even / are unequal 
to zero. This reduces considerably the number of coupled equations to be solved. 
The functions y) satisfy the boundary conditions 



x=0 



and<l>aL(x,?/) 



y=0 



0. 



The partial wave version of the hard-core conditions (|25|) is given by '^aL{,x, y) 
is, by 

^aL{c, y) + y2 dri h^^, (c, y, r]) <^a'Lix\ y') = . 



(38) 
0, that 

(39) 



For the bound-state problem one requires that the functions ^aL{x,y) are square inte- 
grable in the quadrant x > 0, y > 0, i.e., they must satisfy the condition £ -^^2(1^+) 
which follows from (^). A more detailed and useful in bound state calculations is the 
asymptotic condition 



$aL = ^tpiA^) hx{^/E - ei^^y) [aaL,u + o{y ^Z^)] 

V 



+ 



VP 



[AaLie)+o{p- 



(40) 



where E {E < 0) is the bound-state energy, a = {/, A}, and ipi^uix) is the two-body partial 
wave function corresponding to a u-th bound state e; jy for the angular momentum value /. 
Here hx is used for the spherical Hankel function. The coefficients aaL,u and AaL{d) describe 
contributions into $aL (and "^ai) from the (2 + 1) and (1 + 1 + 1) channels respectively. The 
formula (|40|) follows from the asymptotic expression of the total Faddeev component of the 
bound-state wave function (see Ref. f^. Chapter IV, final subsection of §3) which is also 
valid for the hard-core model. 

The asymptotic boundary conditions for the partial Faddeev components of the (2 + 1 — *^ 
2 + 1; 1 + 1 + 1) scattering wave function as X ^ 00 and/or y ^ 00 follow from 
These are 

+ Y.i^,A^)hx'WE-evyy) ^^LyiP) + o [y-"') 

v' 

exp(iA/Ep + i7rL/2) 



(41) 



+ 



Vp 



4'^(P,^) + o(p-^/^) 



where p = |p| is the relative moment conjugate to the Jacobi variable y and the scattering 
energy E is given by = ei^y +p^. The ja' stands for the spherical Bessel function. The 
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value al^iL i,/ represents, at E > evy, the partial amplitude of an elastic scattering, a' = a 

and u' = u, or rearrangement, a' ^ a oi u' ^ z/, process. The functions Al^,'''^{6) provide 
us, at ii^ > 0, the corresponding partial Faddeev breakup amplitudes. Finally the physical 
partial breakup amplitudes are written as 

(^) = 4'^ (^) + E dr^h^a'a'\^.y.r,)A^Ae') (42) 
where 9 and 6' are given by (^). 



IV. APPLICATION TO THE THREE- ATOMIC ^He SYSTEM 



We employed the Faddeev equations (^3|) and the hard-core boundary condition (|39D 
to calculate the binding energies of the Helium atomic trimer and the ultra-low energy 
phase shifts of the Helium atom scattered off the Helium diatomic molecule. As a ^He-^He 
interatomic interactions we use the HFDHE2 |Q and HFD-B [^] potentials of Aziz and 
co-workers. Both HFDHE2 and HFD-B potentials have the form 



y(r)=£|Aexp(-«C + /5r)- 
where C, = r/r^. The function F{Q is given by 

F{C) = 



Cs Cio 



m (43) 



exp[-{D/C-l)Y, ifC<D 
1, ifC> D. 

For completeness the parameters for both HFDHE2 and HFD-B potentials are given in 
Table I. 

In the present work we restrict ourselves to calculations for S-states only. The partial 
components $ao can be obtained in this case from the addition of even partial waves / and 
A with / = A. To demonstrate the feasibility of our formalism and the accuracy which can 
be achieved, we obtained solutions with / = and / = 2 which can be compared with other 
results in the literature. The finite-difference approximation in the polar coordinates p and 
6 has been used for this purpose, a description of which is given in the Appendix. 

Both potentials considered, produce a weakly bound state of the Helium dimer. In our 
calculations we use the value h'^/m = 12.12 K A^. With this value we found the dimer energy 
ea was equal to — 0.8301 mK in the case of the HFDHE2 and to — 1.6854 mK in the case of 
the HFD-B potential. These results are in agreement with other theoretical results found in 



the literature [p9| , p3| , p2| . The estimated experimental value is ~ — ImK [^,^). As to 
the "^He atom-^He atom scattering length, we found that it is 124. 7 A for the HFDHE2 and 
88.6 A for the HFD-B potential. 

Since the Helium dimer bound state exists only in the / = state, the three-body bound 
state boundary conditions (EDI) for the L = channel reads 



^m{x,y) = Sioipd{x) exp{iy/Et - eay) [ao + o (y ^Z^)] 
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where Et and ed are the trimer and dimer energies respectively (expressed in units of A ^) 
and ipd{x) stands for the dimer wave function. 

The results of the Helium trimer ground-state energy Ef'' calculations are presented 
in Table II. It is seen that they are in a good agreement with other results given in the 
literature. Although the two potentials used differ only slightly, they produce important 



differences in the ground-state energy. This is in agreement with the finding of Ref. |32| but 



in disagreement with the statement made in Ref. . It should be further noted that most 
of the contribution to the binding energy stems from the / = A = and / = A = 2 partial 
component the latter being more than 35 %. The contribution from the / = A = 4 partial 
wave was shown in 0] to be of the order of a few per cent. 



In Ref. [^] Cornelius and Glockle investigated the possibility of having Efimov states 
in the Helium trimer. Their work was motivated by the fact that the dimer energy is 
very close to the three-body threshold. Employing the HFDHE2 potential, these authors 
found an excited state at Ef^ = — 1.6mK. This finding was recently confirmed by Esry et 
al. [BH] who also located an excited state at = — 1.517 mK using the same HFDHE2 



inter-atomic interaction. Note that the approaches used in |29[ and |Q are completely 
different. The former is based on the Faddeev integral equations in momentum space while 
in the later on the hyperspherical adiabatic approach. In [^| the improved LM2M2 ^He-^He 



potential was also employed and an excited state at El = — 2.118 mK was found. We 
have also found that the Helium trimer can form an excited state with both the HFDHE2 
and HFD-B potentials. The excited state is present even when only the / = A = partial 



wave is taken into account. This is in agreement with the finding of Ref. |£9|. Our excited 
state results are given in Table HI. By noting that the three-body excited state disappears 
when the interaction strength increases, both, Cornelius and Glockle and Esry et al., identify 
this state as an Efimov one. We have checked a presence of this phenomenon in the case of 
the HFD-B potential and only the partial wave / = A = taken into account. Multiplying 
this potential by an increasing factor (7 > 1 we found that at the beginning the distance 
6 = ed — eI^\ 6 = 6{g), between the trimer and dimer energies E[^\g) and ed{g) increases 
(see Table IV) but thereafter (for A > 1.04) the 5{g) monotonically decreases. As can 
be seen in Table IV, at g ~ 1.19 the value of 5{g) tends to zero, i.e., as in p9| , |33[] , the 
excited state disappears (being covered by the continuous spectrum and probably becoming 
a resonance). This implies the Efimov nature of the excited state energy e[^\ We also 
performed calculations for a Helium atom scattered off a Helium dimer, at L = 0. For this 
we used the asymptotic boundary conditions (^) which, for the L = channel, read 

^m{x, y,p) = dio^pdix) {sm{py) + exp{py) [ao{p) + o (y^^^^)] } 



y P 

The S'-state elastic scattering phase shifts So{p) are then given by 

where So(p) = 1 + 2iao(p) is the (2 + 1^2 + 1) partial component of the scattering matrix. 
Here, we understand the branch of the logarithmic function in such a way that the natural 
"normalization" 6q{0) = 2tt holds. 
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The phase shifts results thus obtained are given in Tables V, VI. We considered incident 
energies below as well as above the breakup threshold, i. e., for the (2 + 1 — > 2 + 1) and the 
(2 + 1 — > 1 + 1 + 1) processes. In order to obtain converged results we were compelled to 
integrate upto a maximum pmax = 400 — 600 A. This comes as no surprise since the two-body 
binding energy is very small implying an extended ^He dimer system and thus the trimer 
wave functions attain their asymptotic values at very large distances. Changing pmax in the 
range 400 — 600 A produces minimal effects on the phase shifts. Such a cut-off radius in the 
three "^He atom problem may be compared with the characteristic values pmax =20 — 30 fm 
for the three-nucleon problem (see, for example, the recent paper and Refs. therein). 
This is not an unexpected result since the dimer wave function generated, for example, by 
the HFD-B potential, behaves as exp(— 0.012 x) at large distances (where x is measured in A) 
while the deuteron wave function as exp(— 0.23a:) (where x is measured in fm). Thus, even 
scaling considerations imply that 20 — 30 fm in the n — d scattering problem are equivalent 
to 400 — 600 A in the three ''He atoms scattering. 

Our estimation for the Helium atom - Helium dimer scattering length 

2 p 

with the HFD-B interactions is 170+5 A in the case if only the / = A = partial Faddeev 
component is taken into account and 145+5 A if the two partial waves with Z = A = and 
/ = A = 2 are considered. In literature we found for 4c only the result 4c = 195 A of Ref. p8 



obtained within a zero-energy scattering calculation based on a separable approximation for 
the HFDHE2 potentials. 

It is interesting to compare the figures for 4c with respective inverse wave numbers 

for the trimer excited state energies. The value of x is given by x = 2\J [ed — e['^'^)/?> where 

both the and have to be measured in A~^. Using data of Table HI for the case of 
the HFD-B interaction we find x~^ 102 A if the only partial wave with / = A = taken 
into account and x~^ 89 A if the two waves with / = A = and / = A = 2 are employed. 
These values are about 1.7 times smaller then the values above for the ^He atom - ^He dimer 
scattering length found on the base of the phase shift results. The situation differs totally 
from the ^He two-atomic scattering problem where the inverse wave number x~^ = 84.8 A 
is rather a good approximation for the ^He-^He scattering length ^^sc = 88.6 A mentioned 
above. Such a significant difference between 4c and x~^ in the case of the ^He three-atomic 
problem may be naturally explained by the Efimov nature of the trimer excited state which 
means that the effective range tq for interaction between ^He atom and '^He dimer is very 
large as compared to the ^He diatomic problem. Unfortunately, insufficient accuracy of 
results for the amplitude ao(p) at p ~ which we have at the moment does not allow us to 
extract the values for the tq. 

We present also a number of figures providing a visual information about the Helium 
atom - Helium dimer scattering. The energy E = +1.4 mK situated above the three-body 
threshold has been chosen for this purpose. In Fig. 2 we plot the Faddeev breakup amplitude 
Aii{e)=Aiio{6) for / = 0. The amplitude ^22(6^) is extremely small, 1^22(6') P < 5-10"^rad"\ 
and therefore it is not shown. The corresponding physical breakup amplitudes Aii{6) = 
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^«zo(^), ^ = 0,2, are plotted in Fig. 3. Both figures explicitly exhibit the importance of the 
inclusion of the / = A = 2 partial waves. The large difference between the results obtained 
for the amplitudes Aix{9) with (/ = A = 0) and (/ = A = 0,2), implies that higher partial 
waves may also be of importance and should be incorporated into the solution in a future 
research. These findings came as no surprise due to the hard-core nature of the underlying 
two-body forces which remain the same for all partial waves. These forces generate strong 
two-body correlations and enhance the role played by higher partial waves. 

As can be also seen in Fig. 3, the breakup is rather suppresed in the vicinity of the 
direction 6 = 30° while the direction of its maximum lies around 6 = 60°. This behaviour 
can be understood from the location of the domains where particles of a particular pair are 
close to each other and thus the elastic scattering of the third particle from the bound pair 
is dominant. 

The absolute value of the Faddeev components ^m{x, y,p) for E = +1.4 mK and / = A = 
is plotted in Figs. 4 and 5 while for / = A = 2 in Fig. 6 and 7 (the corresponding figures 
for the partial wave functions \E'//o(a;, y,p) at the same energy E = +1.4 mK can be found in 
Ref. 1^). In Fig. 4 one can explicitly observe the incident "bound" plane wave ipdix) sm{py) 
which is dominant in ^oqo{x, y,p) at large distances y. However, such a behaviour is not 
present in the partial Faddeev component ^22o{x, y,p) which differs essentially from zero 
only in the vicinity of the triple collision point (see Fig. 6) as in this domain the coupling 
between the channels / = A = and / = A = 2 due to the integral terms of Eqs. (|3^) is 
large. The rather complicated structure of the partial Faddeev components ^iio{x,y,p) in 
the region where the attractions of all three He-He interactions are strong is shown in detail 
in Figs. 5 and 7. It should be noted that the structure of the functions ^iiq{x, y,p) (and also 
of "ifiio^x, y,p) [^) practically does not depend on E in the energy range, ea < E < 2.4 mK, 
considered (compare, for example, our Fig. 5 with Fig. 2 of Ref. p]|). 

V. CONCLUSIONS 

In this work we employ a formalism which is suitable for three-body calculations with 
hard-core potentials. The approach is a hard-core variant of the Boundary-Condition Model 
and, unlike some competing methods, is exact and ideally suited for three-body calculations 
with two-body interactions with a highly repulsive core which can be treated as a hard-core. 
Furthermore the method is feasible not only for bound-state but for scattering processes as 
well. 

We employed the formalism to calculate the binding energy of the '^He-trimer system. 
The results obtained with two realistic ^He-*^He potentials compared favorably with other 
results in the literature. Furthermore, we successfully located an excited state which can 
be interpreted as an Efimov state. This clearly demonstrates the reliability of our method 
in three-body bound state calculations with hard-core potentials. We also endeavored to 
calculate, for the first time, the ultra-low energy scattering phase shifts corresponding to a 
^He atom scattered off a ^He dimer and breakup amplitudes. Using the phase shift results 
we gave as well an estimation for the respective scattering length. 

The effectively hard-core inter-atomic potential together with other characteristics of 
the system, make such calculations extremely tedious and numerically unstable. However, 
this is not the case with our approach where the hard-core is taken from the beginning 
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into account in a mathematically rigorous way. The successful application of the proposed 
method revealed that this method is ideally suited for calculations in systems where the 
strong repulsion in the pairwise forces gives rise to strong numerical inaccuracies which 
make calculations for these molecules cumbersome. Thus the formalism paves the way to 
study various ultra-cold three-atomic systems, and to calculate important quantities such 
as the cross-sections, recombination rates etc. 
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APPENDIX: NUMERICAL METHOD 



The finite-difFerence approximation in polar coordinates p and 9 has been used to solve 
this problem. For this, the grid knots were chosen to be the points of intersection of the arcs 
p — Pi, i — 1,2, ... , Np, and the rays 9 — 9j, j — 1,2, ... , Nq. The value of the parameter 
c ( "diameter" of the particle cores) is chosen in such a way that any further decrease of it 
does not affect the trimer ground-state energy. In the present case a four figure accuracy 
has been achieved with c = 0.7 A. 

The pi points were chosen according to the formulas 



Pi-^ c, 1^1,2,. ..,N^\ 



where Nc''' stands for the number of arcs inside of the core domain and 



The non-hnear monotonously increasing function /(r), < r < 1, satisfying the conditions 
/(O) = and /(I) — 1 was chosen in the form 

in the case of the ground-state calculations and in the form 



f{r) 



aoT , T e [0,ro] 
aiT + T" , T e (to, 1] 



in the case of scattering and excited state calculations. A typical value of the "acceleration" 
a, a > 0, which is satisfactory in ground-state calculations is a = 0.4 (for p^r^ < 100 A). The 
values of ckq, ckq > 0, and a\, a\ > 0, are defined via tq and u from the continuity condition 
for /(r) and its derivative at the point tq. A typical value of tq is 0.2. The value of the 
power u depends on the cut-off radius p^p =200 — 600 A its range being within 3.3 — 4.75. 

The knots for j = 1, 2, . . . , iV^ - were taken according to 9j = arctan(i/j/c). The 
rest knots 9j, j = Np — Nc'^^ + 1, . . . , Ng, were chosen equidistantly. Such a choice of the grid 
is prescribed by the need to have the points of intersection of the arcs p = Pi and the rays 
9 — 9j with the core line x = c as its knots. Furthermore, the grid must be constructed in 
such a manner so that the density of the points is higher where the Faddeev components are 
important, i. e., for small values of p and/or x, and lower in the asymptotic region. Usually 
we took the same numbers of grid points for both 9 and p, Ng = Np. For Nc we chose 
N^'^ = 5. 

The maximal p value used, pmax = PNp, in our ground-state of the Helium ^He trimer 
calculations was 60 A. Beyond this radius the effects on the bound state are minimal. For 
the excited state calculations we were obliged to increase the p^ax to 200 — 400 A and for the 
scattering calculations to 400 — 600 A. 
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A description of a finite-difference algorithm of solving the Faddeev differential equa- 
tions for conventional potentials was given in 0,^. A generalization of this algorithm 
to the boundary-condition model for the three-nucleon problem was previously employed 
in Here, we shall describe in more detail an extension of the algorithm to 



the hard-core boundary conditions problems. For simplicity we restrict ourselves to the 
(2 + 1^2 + 1; 1 + 1 + 1) scattering and the bound-state boundary- value problems where 
only one Faddeev partial equation with / = is considered. 

In the scattering problem, we firstly, in the component y,p) = $000(2;, y,p) explicitly 
separate the initial-state wave function y,p) = ipd{x) sm{py) (see (^5])). As a result, ( |53D 



and ( p9D are reduced to inhomogeneous equations for the remainder = $ — % which differ 
in form from (|33D and (^) only by the presence on the right-hand side of inhomogeneous 
terms F^{x,y) and F'^{y), respectively, whose explicit form is obvious. 

On a fixed arc p = pi oi the polar grid concerned, the values of the function $' and 
inhomogeneous terms F'^(x,y) and F'^(y) form vectors X^'^^ G C^*, JF^*) G M^*, having 



components Xj^^ = ^' {pi cos 6 j, pi sin 6 j) and JF^*-* = F^ {pi cos 6 j, pi sin 6 j) if pi cos 6 j ^ c ox 
Tf' = F'^{pi sin 6j) if pi cos 6j = c. The set of vectors JF^'), i = 1,2, . . . , Np, determines 

the vectors X G C^"" and T G Noo = NoN„: X = ® X'^'l T = © T^'\ In such 
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a representation, Eqs. ( |55D and(^) assumed the form 



A'(o) = 

L,A'(*-i) + (Mi-E/i)A'«+i?iA'(*+i) =J^«, i = \,2,...,Np. ^ ' 

Here, Li, Mi, Ii and Ri are matrices of rank Nq. The matrices and Ri are generated only 
by the radial part of the Laplacian in (p3|) and are therefore diagonal. The non-diagonal 
matrix Mi describes the contribution of the central terms of the radial part of the Laplacian, 
of its spherical part, the potential, and the integral operator on the arc p = pi . In the cases 
where i,j are such that pi cosOj = c, the respective rows of the matrices Li, Mi, and Ri 
are generated by the condition (pOl) . The matrix Ii differs from the unity one only in a row 
corresponding to the boundary condition (|5P|). This row in Ii has zero elements. 

The system ( |Al|) includes Ngp equations for Ngp + Ne unknowns. An additional relation 
that selects a unique solution of ([7n|) follows from the asymptotic conditions 



= B^i^X^""^^ + ao(p)/^,P(^'') (A2) 
where Bjsi^ = diag{6i, &2, • • • , b]\fg} is a diagonal matrix with elements 



bj — ^N, 



(+) 



1 + oipNl ) ' ^No = \l — exp[i/E(p7v,+l - PTVp)] , 

PNp+1 



and V^^p\ V^^p^ G C^», is a vector with components vf"'^ = XiipNp+i,dj) - bjXi{pNp,dj) 
where Xi(P;^) =4'd{p cos 9) exp{ip psinO). 

The condition ( |A2| ) allows the elimination of X^^p'^^^ and reduces the last equation of 
the system to 

L^^A'C^p-i) + (M^^ _ E/^JA^^^") = ^(^"^ + ao(p)^(^'') (A3) 
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where the matrix M^^ and the vector jF^^p) are given by M^^ = Mn^ + R^^B^plNp and 



The system ( |Al|) , after replacing its last equation with (|A3|) , can be written in the form 

{K - Ei)X = + ao(p)J^' (A4) 
where i^' is a three-block-diagonal matrix constructed of the blocks Lj, Mj (or M^^ iii = Np), 



No 



and Ri, i = 1,2, 



,Np, while /, / 



Jj , is a diagonal matrix. Both K and / are 



i = 1 



matrices of rank Ngp. From ( |A1| ), it follows that K has a band structure with band width 



No 



r^'^ with r 



i = 1 



2Ne + 2. The vector P in (^) reads as P 

and^'^^") =^(^''). 

The solution of (^) can be expressed as 

X = Xq + ao(p)A'i 

where the vectors Xq and Xi are determined from 

{K - EJ) A'o = ; (is: - Ei)Xi = T' 



0, i 



1,2,. ..,Np 



(A5) 



(A6) 



in which the inhomogeneous terms are known. 

Having determined the vectors Aq and X\, we can then proceed, via the asymptotics (^5|), 
to find the elastic scattering amplitude ao(p). For this we may use two methods. In the 
first one, we compare the representations (^) and ( [X5| ) on the arc p = Pn in those knots 
{PNp,9j) of the grid for which the value of p^p cos9j belongs to a vicinity of the point Xq 
where the dimer wave function tpd{x) is maximal, ipd{xo) = max?/'(i(a;). In this vicinity, the 
term with the spherical wave exp{i\/E p) / y/p is much smaller than the "elastic" wave term 
ipd{x) exp{ipy) (if pjVp is sufficiently large). Therefore, omitting the breakup term we obtain 
from (|A5|) the following expression 



X, 



(Np) 



xi{Np,ej)- x^ 



{Np) 



(A7) 



where the index j corresponds to the angles 6j for which p^p cosOj ^ xq. 

In the second method we compare the components of ( [A5| ) with the asymptotic rep- 
resentation ( ^5] ) on two successive arcs p = pNp-i and p = pNp, omitting only the terms 
ijjd{x)o{y'^^'^) and exp(iv^p)o(p~^). As a result we find 



with C 



ao(p) 



(-) 

No 



X, 



(Np) 



'-'No 



X, 







(Np) 



Xi{PNp,0,)-C'^-p^\ 4 



ANp-1) 



Xi{p 



Np-l,Oj, 



(A8) 



PNp-l 

PNn 



exp[iV^(pAr — PNp-i)]- As in ([A7|) , the index j corresponds to 



vicinity of the point Xq where the dimer wave function acquires a maximal value. 
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Having calculated ao(p) via ( [X7| ) or (|A8|) we can find, using (|X5|), the vector X^^p^ 



corresponding to the values of the desired function $' on the final arc p = p^p, 
^'{pNp cos 9j,pNp sin6'j) = Xj^''\ and then determine the Faddeev breakup amplitude 

Aoooidj) = Xj'^"^ - ao{p)xiipNp,Oj) y/p^pexp{-i\^pNp) . 

In the bound-state problem we deal with the same system of equations ([Al|) for Xj^^ = 
^{pi cos dj, Pi sin 9j) where now y) stands for a bound-state wave function satisfying the 
asymptotic conditions (0). Of course the inhomogeneous terms JF^*) vanish in this case. 

To eliminate the vector X'^^p^^^ from the last (z = Np) equation of ( |AT|) we use the 
representation P^). For angles corresponding to the knots of the arc p = PNp lying inside 
the core domain, p^^cosdj < c, we write the components Xj^''^ and X^^'''^^^ on the two 
successive arcs p = p^p and p = PNp+i, taking into account the condition ipd{x) = 0, 
X < c, and neglecting the terms exp{i^/E p)o{p^^) . Then we find Xj^'''^^^ = C^^^X^^''^ . 
For angles 9j corresponding to knots of the arc p = PNp lying outside the core domain, 
PNp cos 6j > c, we write the expression for the components X^^^ on three successive 
arcs p = Pi, i = Np — 1, Np, Np + 1, neglecting the terms 4'd{x)o{y~^) and exp(i\/Ep)o(p~^). 
Using this expression for i = Np ~ 1 and i = Np we can express ao and A{6j) in terms of 
Xj^'' ^'^ and Xj^''\ Then we may express Xj^'''^^'^ in terms of X^^'' and Xj^''^ using (^if ) 



for 2 = iVp -|- 1. Thus, finally, the last equation of (|A1|) becomes 



V^"""-'^ + {Mnp - EiNp)x^''p'> = 

where the matrices L^^ and M^^ are given by 

Lnp = Lnp + RNplNpWNp , Mnp = Mnp + i?7v, /tv, W^TV, • (A9) 

The W = dia:g{w[^''\ . . . ,w\!^J'^} and W = diag{w[^''\ . . . ,w^j^J'^} are diagonal matrices 
with 



Xi{PNp+i, Oj ) - C^Np X^ (PNp ' ) 
(Np) I r-3i , PNpCOSej>C, 

^ =) xi{PNp^i,o,)-ci^ xi{PNp,e,) (Alo) 



, PNp COS 9j <c, 



and 



I Cl^pXiipNp-i^dj) - Cl'' Xl{pNp+l,Oj) 

J n^i ' pNpCose,>c, 

~ xi{PNp-u0,)-c^Np^ xi{PNp,e,) (All) 

Cn„^ , pNp COS 9 j < c. 



where now Xi(P) ^) = i^dip cos 9) exp(i V E — edP sin 9) . Note that the matrices Lnp and Mnp 
depend on the energy E since the function xi and the coefficients C^^J are functions of E. 
Therefore the total matrix K of the system obtained is also a function of E, K = K{E). 
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In this work we searched for binding energies of the ^He trimer as roots of the determinant 
det{K{E) - EI}. 

The use of the asymptotic boundary conditions ( PI) in the form of ([A9|)-( |70A| ) is ex- 
tremely important when searching for the excited Ej:^^ state. It is difficult to locate this state 
if the term a.oijjd{x) exp{i\/ E — e^y) is omitted. This means that the dimer wave function 
ipd{x) gives a decisive contribution into the excited state. In contrast, omitting this term 
(and replacing wf''\ wf"'^ given by ( |ATO|) and (|ATID with wf"'^ = 0, wf'^ = C^^J for 
all j = 1, . . . ,Ng) in the ground-state calculations simplifies the problem considerably by 
allowing the decrease of the cut-off radius to 60 A. Otherwise, to obtain the correct result 
for E^^^ we had to increase pn up to 150 — 200 A as the dimer wave function falls off slowly. 
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TABLES 

TABLE L The parameters for the ^He— ^He potentials used. 



Parameter 



HFDHE2 g 



HFD-B H 



e(K) 
(A) 
A 
a 
P 
Ce 
Cs 
Cio 
D 



10.8 
2.9673 
544850.4 
13.353384 


1.3732412 
0.4253785 
0.178100 
1.241314 



10.948 

2.963 

184431.01 

10.43329537 

-2.27965105 

1.36745214 

0.42123807 

0.17473318 

1.4826 



TABLE IL Bound state energy results for the Helium trimer. The (absolute) values of 
E^^^ are given in K. The grid parameters used were: Nq = Np = 275, a = 0.4, and /Omax =60 A. 



Potential 


Faddeev equations 


Variational 
methods 


Adiabatic 
approach 


/ This work |3Cl |29| |28| 


[22| |2C| 




HFDHE2 


0.084 0.082 0.092 


0.1173 


0.098 


0,2 0.114 0.107 0.11 


HFD-B 


0.096 0.096 


0.1193 




0,2 0.131 0.130 



TABLE III. Excited state energy E^ results for the Helium trimer. The (absolute) values of 
E^-^^ are given in mK. The grid parameters used were: Nq = Np = 252, tq = 0.2, u = 3.6 and 

Pmax =250 A. 



Potential 


I 


This work 


|29| 


|28| 


11 


HFDHE2 





1.5 


1.46 


1.04 


1.517 




0,2 


1.7 


1.6 






HFD-B 





2.5 










0,2 


2.8 
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TABLE IV. Dependence of the dimer bound state, €4, and trimer excited state, E^. , energies 
on the multiplying factor (potential strength) g in the case of the HFD-B ^He-'^He interaction. The 
values of energies e^, E^^"^ and their difference 5 = — E^^^ are given in mK. The grid parameters 
used were: Ng = Np = 550, tq = 0.215, u = 4.75 and pmax =350 A. 



9 


ed 


El" 


6 


0.975 


-0.036 


-0.308 


0.272 


1.00 


-1.685 


-2.485 


0.800 


1.04 


-9.368 


-10.353 


0.985 


1.10 


-32.222 


-32.777 


0.556 


1.16 


-68.150 


-68.334 


0.184 



TABLE V. Phase shift results for the L = 0, l = X = partial wave obtained with the HFD-B 
^He-^He potential. The grid parameters used were: Nq = Np = 320, tq = 0.2, v = 4.5, and 

=400 A . 



E (mK) 


5q (deg) 


E (mK) 


do (dcg) 


E (mK) 


5o (dcg) 


-1.68535 


359.2 


-1.4 


315.8 


0.1 


276.8 


-1.6853 


358.9 


-1.25 


309.0 


0.4 


272.6 


-1.685 


357.5 


-1.1 


303.6 


0.7 


268.9 


-1.68 


352.3 


-0.95 


299.1 


1.0 


265.7 


-1.67 


347.2 


-0.8 


295.0 


1.4 


261.8 


-1.66 


344.2 


-0.7 


292.5 


1.6 


260.0 


-1.65 


341.4 


-0.5 


287.9 


1.8 


258.4 


-1.60 


333.4 


-0.4 


285.7 


2.1 


256.2 


-1.5 


322.3 


-0.2 


281.7 


2.4 


254.2 



TABLE VI. As in Table V but with the inclusion of the I = \ = 2 partial wave. 



E (mK) 


5q (deg) 


E (mK) 


5q (deg) 


E (mK) 


6q (deg) 


-1.68535 


359.3 


-1.4 


323.0 


0.3 


283.7 


-1.6853 


359.0 


-1.1 


312.6 


0.7 


278.7 


-1.685 


357.8 


-0.8 


304.6 


1.0 


275.4 


-1.68 


353.3 


-0.6 


299.8 


1.4 


271.5 


-1.60 


336.4 


-0.4 


295.5 


1.8 


268.0 


-1.5 


328.0 


-0.1 


289.7 


2.4 


263.5 
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FIGURES 




A I ^. (rad ^ ) 



0.010 - 



0.005 - 




0.000 ' ' ' ' ' ' ' ' ^"^ 1 ' ' ' ' ' ' 

10 20 30 40 50 60 70 80 90 

(degree) 

FIG. 2. The square of the modulus of the Faddeev breakup amphtude Aoo{6) for HFD-B 
^He-^He potential at £^ = +1.4 mK. Curve 1 corresponds to the L = 0, l = X = partial wave 
while curve 2 was obtained with the inclusion of the L = 0,l = X = 2 channel. The grid parameters 
used were the same as in Table V. 



0.015 - 



0.010 - 




(degree) 

FIG. 3. The squares of the moduli of the physical breakup amplitudes v4oo(^) (curves 1, 2) 
and A22{0) (curve 3) for the HFD-B ^He-^He potential at £^ = +1.4 mK. Curve 1 corresponds to 
the inclusion of the L = 0, Z = A = channel only, while curves 2 and 3 were obtained with the 
inclusion of both Z = A = and Z = A = 2 partial waves. The grid parameters used were the same 
as in Table V. 
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0.119 

0.0952 

0.0714 

0.0476 




FIG. 4. Absolute value of the Faddeev component ^ooo{x,y,p) for the HFD-B ^He-^He poten- 
tial at E = +1.4 mK. The grid parameters used were the same as in Table V. The values of x and 
y are in A. 



0.0593 

0.0475 

0.0356 




15 10 



X 

FIG. 5. Detail of the \^ooo{x, y,p)\ surface shown in Fig. 4. 
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0.00673 
0.00538 
0.00404 
0.00269 
0.00135 




FIG. 6. Absolute value of the Faddeev component ^22o{x, y,p) for the HFD-B ^He-^He poten- 
tial, at £^ = +1.4 mK. The grid parameters used were the same as in Table V. The values of x and 
y are in A. 




FIG. 7. Detail of the \^22o{x, y,p)\ surface shown in Fig. 6. 
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